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Abstract The recent discovery of high frequency os- 
cillations during giant flares from the Soft Gamma Re- 
peaters SGR 1806-20 and SGR 1900+14 may be the first 
direct detection of vibrations in a neutron star crust. If 
this interpretation is correct it offers a novel means of 
testing the neutron star equation of state, crustal break- 
ing strain, and magnetic field configuration. We review 
the observational data on the magnetar oscillations, in- 
cluding new timing analysis of the SGR 1806-20 giant 
flare using data from the Ramaty High Energy Solar 
Spectroscopic Imager and the Rossi X-ray Timing Ex- 
plorer. We discuss the implications for the study of neu- 
tron star structure and crust thickness, and outline areas 
for future investigation. 
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■ 1 Introduction 

The Soft Gamma Repeaters (SGRs) are thought to be 
magne tars, neutron stars with magnetic fields in excess of 
10 14 G ijDuncan fc ThomDSOiJl992lThonroson fc Duncanl 
Il995|) . Decay of the strong field powers regular gamma- 
ray flaring activity that culminates, on rare occasions, in 
a giant flare with a peak luminosity in the range 10 44 - 
10 46 erg s _1 . The three giant flares detected to date con- 
sist of a short, spectrally hard initial peak, followed by 
a softer decaying tail that lasts for several hundred sec- 
onds. Pulsations with periods of a few seconds are visible 
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in the tail and reveal the neutron star spin. Their pres- 
ence is thought to be due to a fireball of ejected plasma, 
trapp ed near the stellar surface b y the strong magnetic 
field ijThompson fe Duncadll995|) . 

Powering the giant flares requires a catastrophic global 
reconfiguration of the magnetic field. The coupling be- 
tween the field and the charged particles in the neutron 
star crust means that this is likely to be associated with 
large-scale crust fracturing ^Flowers fc Rudermanlll977t 
iThompson fc DuncaiJll993. 120011: ISchwartz et al.ll2005[) 
This in turn should excite global seismic vibrations: on 
Earth seismologists regularly observe global modes af - 
ter large earthquakes (see for example lPark et al] l2005). 
In the SGR case, the coupling of field and crust should 
cause the modes to modulate the X-ray lightcurve. 

Various different types of oscillation are possible, but 
theory suggests that the easiest to excite an d observe 
shoul d be the toroidal shear modes of the crust ijBlaes et al.1 
1989). The particular harmonics excited will depend on 
the size, shape and speed of the fracture. The observed 
mode frequencies depend on the neutron star mass and 
radius (via gravitational redshift and the influence of 
the equation of state on crust structure), crustal com- 

position, and magnetic field strength and configuration 

( Hans^n^^CjoffiLgj^lSchu^ 

1988; Strohmayer 1991; Duncan 1998; Messios, PapadoDOulous & Stcrgioi 



200lUPird i2005). Detection and identification of crustal 
modes would therefore probe all of these areas of neutron 
star physics. 

In 2004, SGR 1806-20 emitted the m ost powerful 
flare e ver recorded. In a landmark paper, llsrael et alJ 
(2005) reported the detection of Quasi-Periodic Oscilla- 
tions (QPOs) in the 18-93 Hz range in the decaying tail of 
the fla re. Prompted by this result. Istrohmaver fc Watts! 
(|2005l) re-analysed the 1998 giant flare from SGR 1900+14 
and found QPOs in the range 28-155 Hz. In both cases 
tentative identifications can be made with a sequence 
of toroidal modes. More in-depth analysis of the SGR 
1806-20 flare has since revealed additional QPOs, in- 
cluding the first evidence for a higher radial overton e 
(|Watts fc Strohmaveill2006t IStrohmaver fc Wattsl 2006) . 
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We review all of the observational results in Section El 
In Section we discuss the results in the context of the 
toroidal mode models, and show how they can be used 
to constrain neutron star parameters including the crust 
thickness. We conclude in Sections 0] and El with a dis- 
cussion of outstanding issues. 



2 Observational results 

2.1 SGR 1806-20 

The December 2004 flare was the brightest ever recorded, 
with a peak luminosity of ~ 10 46 erg s _1 . Analysis of 
data from the Ros si X-ray Timing Explorer (RXTE) by 
llsrael et "all l|2005|) revealed a transient 92.5 Hz QPO in 
the tail of the flare, associated with a particular rota- 
tional phase. The QPO appeared around 170s after the 
main peak, at the same time as an apparent boost in 
impulsed emission. The presence of weaker 18 and 30 Hz 
features at late times was also suggested. 

The toroidal modes are labelled by the standard quan- 
tum numbers I, m and n, the first two being angular 
quantum numbers, the final one denoting the number 
of radial nodes. T he 30 and 92.5 Hz features found by 
llsrael et all l|2005l) can be identified as the fundamental 
I = 2, n = mode and the I = 7, n = harmonic (see 
Section |21 for a discussion of current mode models). The 
18 Hz feature, by contrast, is at too low a frequency to 
be identified as a toroidal mode of the crust. 

The other spacecraft with high time resolution data 
of the SGR 1806-20 giant flare was the Ramaty High 
Energy Solar Spectroscopic Imager (RHESSI), a solar- 
pointing satellite that covers a wider energy band than 
RXTE. RHESSI's detectors are split into front and rear 
segments. Strong albedo flux from the Earth affected the 
time resolution of the rear segments (making them un- 
suitable for studying frequencies > 50 Hz), but when 
the rear segments are included the countrate recorded 
by RHESSI exceeded that of R XTE. 

IWatts fc Strohmaverll|2006|) analysed the RHESSI data 
of the flare and confirmed the presence of the 92.5 Hz 
QPO, at the sam e time and rotational phase found by 
llsrael et "all l|2005|) . At low frequencies, where the larger 
RHESSI countrate gave added sensitivity, broad QPOs 
at 18 and 26 Hz were found 50-200 s after the main flare, 
at the same rotational phase as the 92.5 Hz QPO. Al- 
though there was a weaker feature at 30 Hz, it was not 
at the 3cr level after accounting for the number of trials 
so we were not abl e to make a robust confirmation of the 
llsrael et "all l|2005|) result. Subsequent closer comparison 
of the two datasets, however, indicates that RHESSI's 
countrate is not the only factor affecting its sensitivity. 
As such we now believe that the 30 Hz detection is robust 
(see below). 



Also detected in the RHESSI data was a 625 Hz 
QPO 1 at an energy band nominally higher than that 
recorded by RXTE 2 . A full discussion of the QPO prop- 
erties and the s ignificance of the detection c an be found 
in Section 2 of IWatts fc Strohmaverl i(2006t) . Compared 
with the 92.5 Hz QPO the 625 Hz QPO is in a higher 
energy band, emerges earlier in the tail of the flare, and 
is at a different rotational phase. Excitingly, however, it 
is also at the approximate frequency expected for the 
radial overtone n = 1 toroidal modes ljPiroll2005|) . 

We have recently completed a more in-depth analysis 
of th e now public RXTE dataset ijStrohmaver fe Watts! 
2006) . We start by discussing the low frequency QPOs. In 
this regime the RXTE countrates are lower than RHESSI's. 
However, at times when the 18 and 26 Hz QPOs are 
particularly strong in the RHESSI dataset they are also 
detected in the RXTE data. We are thus able to con- 
firm both the frequencies and the rotational phase depen- 
dence. Further ana lysis of the 30 Hz QPO confirms the 
detection made bv llsrael et alJ l|2005fl . and reveals that 
this QPO is also strongly rotational phase dependent. 
Rather surprisingly, given the difference in countrates, 
this feature is far stronger in the RXTE data than in the 
RHESSI dataset. There are several factors, however, that 
^could account iOT^this discr epancy (see the discussion in 

Looking to higher frequencies, a rotational phase de- 
pendent search reveals additional QPOs at 150 Hz, 1840 
Hz and, most excitingly, at 625 Hz (the latter appearing 
at the same rotational phase as the 92.5 Hz QPO). Full 
details of the QPO properties and the significance of the 
detect ion can be found in Section 2.1 o flStrohmaver fc Wattsl 
(2006]). There are, it turns out, intriguing differences be- 
tween the RHESSI and RXTE 625 Hz QPOs. The RXTE 
QPO has lower fractional amplitude, lower coherence, 
a different rotational phase association, and appears at 
later times in a lower energy band. One possibility is 
that we are seeing time evolution of one QPO whose 
photon energy and amplitude are falling over the course 

1 Grid shadowing in RHESSI can give rise to spurious high 
frequency signals, and there is a known artifact at 718 Hz 
in the giant flare data (Gordon Hurford, private communi- 
cation). The phenomenon affects certain detectors far more 
strongly than others since it is due to spacecraft geometry. 
It is also (perhaps unsurprisingly) more prominent at lower 
energies. The fact that the 625 Hz QPO is dependent on the 
rotational phase of the magnetar, despite there being com- 
parable numbers of photons at other rotational phases, is a 
good indication that it is not a detector artifact. In addition it 
is only detected at higher energies, where shadowing should 
be less of an issue. However, to be sure we also tested the 
variation of QPO power when photons from individual de- 
tectors were removed from the analysis. The drops in power 
were consistent with the drop in countrate, confirming that 
the signal is not anomalously strong in one or two detectors. 
As such we are confident that the signal is indeed associated 
with the source. 

2 Direct comparison of recorded energies is difficult because 
for both spacecraft the flare was off-axis, resulting in substan- 
tial scattering within the body of the satellite. 
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of the tail. The second possibility is that we a r e see- 
ing two different modes. As pointed out bv iPirol (|2005[) 
the n = 1 modes are nearly degenerate in I, so there 
are several different modes with very similar frequencies. 
Despite the differences, it seems extremely unlikely that 
two independent instruments would detect signals at a 
consistent frequency unless that frequency were intrinsic 
to the source. 



2.2 SGR 1900+14 

In August 1998 SGR 1900+14 emitted a giant flare with 
a peak luminosity ~ 10 44 erg s _1 . The flare was detected 
by RXTE, albeit with data gaps due to the configura- 
tion o f the spacecraft at that time. lstrohmaver fc Wattsi 
l)2005|) started by analysing each good interval separately, 
and discovered a strong transient 84 Hz QPO during a 
~ 1 s interval about 60 s after the main flare. Folding 
up data from more intervals revealed a pair of persistent 
QPOs at 53.5 and 155.5 Hz, with a weaker feature at 28 
Hz, all at the same rotational phase as the 84 Hz signal. 
No QPOs were detected at other rotational phases. The 
scaling of the QPO frequencies is what the existing mod- 
els suggest for the n=0, 1=2, 4, 7 and 13 toroidal modes 
of the crust. 



3 Constraining neutron star properties 

Detection of signals at similar frequencies in the giant 
flares from two different SGRs implies that the same pro- 
cess is operating in both objects. In addition, the strong 
rotational phase dependence of all of the detected QPOs 
provides strong evidence that the modulations are as- 
sociated with the stellar surface. In this section we dis- 
cuss the implications of the QPOs in the light of the 
toroidal mode model. Existing models do have short- 
comings, which we discuss in more detail in Section 
However, crust mode models remain at present the most 
promising mechanism and we will proce ed accord ingly. 

For a non-rotating, nonmagnetic star. lDuncanl {l998) 
estimates the frequency v of the fundamental I = 2,n = 
toroidal mode (denoted 2*0) to be 



^(2*0) 



29.8 t 171 ; - 71 ^ 1 ! 2 Hz 



0.87i?i + 0.13Mi. 4 i? 



(1) 
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where R w = R/W km and M 1A = M/1.4M Q . The fre- 
quencies of the higher order n = modes are given by 



v{ito) = ^(2*0) 



1(1 + 1) 



G 



1/2 



1 



B 



1/2 



(2) 



where the final factor is a magnetic correction, w 
4 x 10 15 p?4 4 G and pi 4 ~ 1 is the density in the deep 
crust in units of 10 14 g cm~ 3 . In deriving this equa- 
tion it is assumed that magnetic tension boosts the field 




1.4 1.6 

Mass (solar masses) 

Fig. 1 Mass and magnetic field required to give the 
toroidal mode frequencies inferred for SGR 1806-20 and 
SGR 1900+14. We s how r esults for four of the EOS from 
ILattimer fc Prakashl 120011k in order of increasing stiffness 
they are WFF1, AP4, AP3 and MS0. The upper line for 
each EOS indicates the parameters necessary to give the SGR 
1806-20 frequencies; the lower line those for the lower SGR 
1900+14 frequency. The horizontal line indicates the uncer- 
tainty in the position of the footpoints of the lines due to the 
width of the QPOs. 



isotropically. Field configuration and corresponding non- 
isotropic e ffects could however alter this correction dra - 
matically l)Messios. Panado pojflous^ fe Stereioulasl 2001'l. 
More recent calculations bv lPirol l)2005T) using better mod- 
els of the neutron star crust confirm these frequency es- 
timates. 

The detection of a set of modes from both SGR 1806- 
20 and SGR 1900+14 with the expected [;(2 + l)] 1/2 scal- 
ing in frequency is indicative of the presence of toroidal 
modes. We can then ask what can be learnt from the fact 
that the fundamental 2^0 mode frequency seems to differ 
for the two SGRs. For SGR 1806-20 it is inferred to be 30 
Hz, whereas for SGR 1900+14 it is lower, at 28 Hz. For 
this to be the case, the properties of the two stars must 
differ. Given an equation of state (EOS), equations (QJ 
|5J| specify the relationship between mass and magnetic 
field necessary to gi ve modes at the inferred freq uencies. 
Figure^ taken from lStrohmaver fc Watts! l)2005|) . shows 
for several EOS the stellar parameters that give 2^0 os- 
cillations at the frequencies inferred for the two stars. 

Several conclusions follow from Figure ^ Firstly, un- 
less the masses differ substantially, SGR 1806-20 must 
have a stronger magnetic field than SGR 1900+14. This 
acc ords with estimates of field strength from timing stud- 
ies l|Woods et al.l2002|) . Secondly, the highest masses and 
the hardest EOS require fields far higher than those in- 
ferred from timing, whereas the softest EOS require high 
masses. Clearly these inferences are all model-dependent, 
and could change as the models get more sophisicated. 
However, it illustrates the potential of neutron star seis- 
mology to constrain the equation of state. 
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If the 625 Hz QPOs detected in the SGR 1806-20 
flare are indeed n — 1 toroidal crust modes then the 
implications are profound, since from the ratio of fre- 
quencies of the n — and n = 1 modes one can deduce 
the thickness of the crust, AR (gravitational redshift fac- 
tors scale out). In the limit of a thin crust {AR <C R) 
and constant shear wave speed llHansen fc Cioffil 1980t 



iMcDermott. van Horn fc Hansen! 1 19881: iPird l2005fl one 

can show that 



AR 
~R 



1/2 



v(it x ) 



(3) 



So given vfaa) = 30 Hz and v{it\) = 625 Hz we find 
AR/R = 0.06. A more sophisticated estimate that does 
not assume a thin crust or a constant shear speed gives 
AR/R in the range 0.1 - 0.12 l|Strohmaver fc WattsfcOOfil) . 

Measuring the crust thickness also gives an additional 
constraint on the EOS (independent of redshift effects), 
since stellar models of a given mass computed with differ- 
ent EOS will in general have crusts with different depths. 
The article by Lattimer in this volume outlines how this 
could be used to constrain the nuclear symmetry energy 
and the nuclear force model. 



4 Theoretical issues 

Neutron star seismology, revealed during giant flares from 
magnetars, has great potential as a probe of stellar struc- 
ture, composition and magnetic field geometry. There are 
many theoretical issues, however, that remain to be re- 
solved. One of the main issues is that of the coupling 
between the crust and the fluid core due to the pres- 
ence of the strong magnetic field. Global magneto-elastic 
modes of the whole star may be necessary to explain, in 
particular, the QPOs at 18 and 26 Hz, which do not sit 
comfortably within existing toroidal mode models. The 
coupling between the crust and the field will also deter- 
mine how the modes modulate the lightcurve. 

Modes of the neutron star crust (most likely cou- 
pled to the core via the magnetic field) remain the most 
promising mechanism identified to date for the QPOs. 
Although several other suggestions have been made, all 
have serio us difficulties-M odes of the magnetosphere, men- 
tioned by iLevinl l)2006f) , are likely to have too high a 
frequency due to the high Alfven speed in this region. 
Modes of the trapped fireball are unlikely since one would 
expect the frequency to change as the fireball shrinks: no 
such correlation is observed. The third possibility, raised 
by Alpar at this meeting, is interaction with a debris 
disk and a mechanism similar to that thought to gen- 
erate QPOs in the accreting systems. However even if 
a small disk were present in these systems, it would be 
difficult to explain the rotational phase dependence. 

The complex temporal variation of the QPOs also 
requires explanation. In both giant flares some of the 



QPOs seem to be highly transient, whereas others per- 
sist for most of the flare. We mentioned in Section 
the possible evolution of the 625 Hz QPO seen during 
the SGR 1806-20 flare. However the 92.5 H z QPO also 
shows variation in frequency and amplitude l|lsrael et alJ 
l2005UStrohmaver fc WaHsl [2006h Evolution of the mag- 
netic field in the aftermath of the flare, relaxation of 
the deep crust, or evolution of the surrounding plasma 
are all candidates for causing such variations. Damping 
and excitation mechanisms are a critical area for future 
study. 



5 Conclusions 

The discovery of high frequency oscillations during gi- 
ant flares from the Soft Gamma Repeaters SGR 1806-20 
and SGR 1900+14 may be the first direct detection of 
vibrations in a neutron star crust. If this interpretation 
is correct it offers a new means of testing the neutron 
star equation of state, crustal breaking strain, and mag- 
netic field configuration. In particular, if the mode inter- 
pretation is secure, it allows us to make the first direct 
estimate of the thickness of a neutron star crust. This is 
particularly impressive when one considers the fact that 
all of the data acquired so far have been obtained by 
chance using satellites that were observing other objects. 
If a rapid-slew instrument such as SWIFT could be con- 
figured to record high time resolution data throughout 
the tails of giant flares from known SGRs, the potential 
for additional discoveries is immense. 
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